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HISTORICAL NOTE ON HEAVISIDE’S OPERATIONAL 
METHOD. 


By N. W. McLACHLAN. 


OPERATIONAL methods in pure mathematics, e.g. the work of Boole, 
Graves and Murphy, existed many years before Oliver Heaviside 
introduced new ideas into mathematical physics. His first paper 
using such ideas was published in the Proceedings of the Royal Society 
in 1893.* Remarks therein indicate his unfamiliarity with some of 
the literature of the subject, owing to lack of reference facilities, 
for example, Liouville’s work on fractional differentiation. Heavi- 
side’s main contributions to operational method seem to be: (a) the 
substitution p= (d/dt)”; (b) expansion in inverse powers of p and 

“ algebrisation ” using the formula p»=t/r(1 +v); (c) incorpor- 
ation of the initial conditions prior to solution ; (d) the expansion 
theorem ; and (e) expansion in ascending fractional powers of p to 
obtain an asymptotic formula. His main purpose was to solve 
transient problems in cable telegraphy and telephony, so he did not 
justify the analytical processes used. His method always yielded (in 
his hands!) the correct result, and he argued that this was in itself 
ample justification for the operational method. He wrote “ The 
best proof is to go and do it”’.+ Moreover, contemporary pure 
mathematicians looked upon his work with scepticism until such 
time as it was put ona sound foundation by Bromwich f and Wagner § 
working independently. The boycotting of his work elicited some 
apposite remarks from Heaviside, as, for example: “ Even Cam- 
bridge mathematicians deserve justice. . . . As regards their want of 
sympathy with less conventional men, it is not sympathy that is 
particularly wanted. ... What one has a right to expect, however, 
is a fair field, and that the want of sympathy should be kept in a 
neutral state, so as not to lead to unnecessary obstruction.” 

“ Orthodox mathematicians, when they cannot find the solution 
of a problem in a plain algebraical form, are apt to take refuge in 
a definite integral, and call it a solution. ... But it may be just as 
hard, or harder, to interpret than the differential equation of the 
problem in question.” || 

This psycho-polemical phase of the subject has passed away (we 
hope!), and Heaviside’s operational calculus is now adopted (some- 
times with considerable modification as will be shown later) as a 
standard technique in the solution of certain classes of problems, 
frequently those involving transient effects and wave motion. In 
the memorable paper where Bromwich ** established the validity of 


* Proc. Roy. Soc. A. 52, 504 (1893). 
+ E. M. Theory, vol. 2, 34 (1922). t Proc. L.M.S. (2), 15, 412 (1916). 
§ Archiv fiir Elektrotechnik 4, Band 5, 159 (1916). 


|| H.M.7T. vol. 2, pp. 10, 11 (1922 edition). We disagree with Heaviside: see 
for example Math. Gaz. 22, 37, 1938, where a compact solution of a cable problem 
includes a definite integral which can be interpreted physically. 


** loc. cit. 


256 THE MATHEMATICAL GAZETTE 


Heaviside’s expansion theorem, he writes, “ The investigations lead- 
ing up to the formula (the expansion theorem) * are most instructive 
and will repay careful study: but I am not sure that they have 
been sufficiently appreciated in the past. It is almost certain that 
few readers have fully grasped the complete and general character 
of the solution : and it is for this reason that I wish to call attention 
to it here.” 

Some forty years after publication of the expansion theorem, one 
of the leading authorities of the Cambridge school of mathematicians 
wrote: “‘ We should now place the operational calculus with Poin- 
caré’s discovery of automorphic functions and Ricci’s discovery of 
the tensor calculus, as the three most important mathematical 
advances of the last quarter of the nineteenth century. Applications, 
extensions, and justifications of it constitute a considerable part of 
the mathematical activity of to-day.” f 

Bromwich ł established the validity of Heaviside’s method of 
operators—at least in so far as it involved the evaluation of residues 
at poles—by complex integration. Amongst other things, he pointed 
out that when f(t) is zero outside the range 0 <<t<7, if 
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where &(p) is to be regarded as a convenient transform for f(t). 
Suppose we solve the equation 


dy dy ; 
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subject to the conditions y’=y=0 at t=0, by the Bromwich 
method. Then we assume that § 
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which on substitution in the left-hand side of (3) yields 
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Since sin é continues indefinitely, r=% in (1), so 


E(p) =f, | ae sin tdt =fy[(p2+1), oasen (5) 


* Electrical Papers, vol. 2, p. 259. 

t E. T. Whittaker, “ Oliver Heaviside ”, Bull. Math. Soc. Calcutta, 20, 199 
(1928-29). 

$ loc. cit. 

§ All the singularities of the integrand lie to the left of the contour. If they 
are all poles, c + ¿œ may be replaced by a circle enclosing them. 
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and by (2) 
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Hence by (4) and (6) 
x(p) =fo/(p? tap +b) (pP? +1), ..ccccrccccccscseecenees (7) 
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which is evaluated in the usual way by the calculus of residues. 

The origin of formulae (1) and (2) is sometimes misquoted by 
writers on the subject of operators, so we give Bromwich’s remarks * 
thereon verbatim : “ The formulae may be deduced, as a matter of 
simple transformation, from the ordinary form of Fourier’s theorem. 
Formulae substantially equivalent were obtained in this way by 
Riemann in his paper on the distribution of primes (Ges. Werke, 
p. 140); and the actual formulae were deduced similarly by Mac- 
donald (Proc. L.M .S. 35, 428 ; 1902). From the point of view of the 
theory of functions of a complex variable, more complete discussions 
have been given by Pincherle and Mellin ; see for instance a paper 
by the latter in Math. Annalen, Band 68, p. 305, where references 
to earlier investigations will be found.” 

Heaviside did not seem to be favourably impressed by the method 
of contour integration, as will be realised on reading the following 
excerpt from a letter t he wrote to Bromwich on 7th April, 1919, 
from Torquay, where he lived till his demise in 1925 : 


“DEAR Dr. BROMWICH, 


Yours 5:4:19, Caesar and Pompey: especially Pompey. 
. What a time it takes! I rejoice to know that you have seen 
the simplicity and advantages of my way. Only a war could have 
done it. Did you see in the paper that I was the cause of the war, 
and that I brought it on to wake up England and especially Cam- 
bridge? Now let the wooden headed rigorists go hang, and stick 
to differential operators and leave out the rigorous footnotes. It is 
easy enough if you don’t stop to worry. As I said and Lord R.{ 
repeated, logic is the very last thing. I never could stomach your 
complex integral method, though I understand your affection for 
it... . I don’t care much about the expansion theorem myself now, 
it is so tame... .” 
After this brief résumé of Bromwich’s work, we turn to Wagner, § 
who published a paper almost contemporaneous with that of the 
former writer. He commenced with the integral 


_ Jo i n2 f t <0 
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* loc. cit. p. 412. 


t The original was given to A. T. Starr when at Cambridge under Bromwich. 
The author is indebted to the former for a lantern slide of the original script. 


+ Presumably the late Lord Rayleigh. § loc. cit. 
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taking the imaginary axis indented on the right at the origin as 
contour. Integral (9) was used to define Heaviside’s unit function. 
Wagner went on to state that Heaviside’s first expansion theorem 
was equivalent to evaluation of an integral of the form 


(30) 
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by finding the residues at the poles of the integrand. No formal 
proof was given, however, but a number of examples from engineer- 
ing practice were treated in a way comparable with that of Bromwich. 

The subject of branch points was not mentioned in either of these 
papers. In fact it was not until some years later that the topics of 
branch points, equivalent contours (e.g. —œ, (0+), —«) and the 
asymptotic expansions used by Heaviside were considered with 
respect to the complex variable. 

The contour used by Bromwich, namely, c-1” to c+i1 is 
more general than that of Wagner, and it is now employed for the 
interpretation of an operational expression ¢(p), as in the following 
integral : 

7 (t) B Iri 
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all singularities of the integrand being to the left of the contour. 
In its relation to operational] procedure, this formula can be con- 
sidered to be based upon (2) and (10). Accordingly it may be 
regarded as the Bromwich-Wagner integral. In modern practical 
applications it is useful for interpreting operational forms ¢(p), 
where 


$(p) =p i TOT cescsssesereeseeseneenes (12) 


When certain conditions * are fulfilled, (11) and (12) constitute a 
particular case of the Mellin inversion theorem.t ¢(p)/p is the 
Laplace transform of f(t). 

The contour c+1œ, where cœ Q is not necessarily constant, 
may for present purposes be regarded as the Bromwich contour. 
Provided the integrals along the arcs in the second and third quad- 
rants, joining the extremities of the contours, vanish as their radii 
—>æ, any contour starting in the second quadrant, crossing the 
real axis at z=c, > 0 and terminating in the third quadrant, is 
equivalent to c+to. In certain cases the equivalent contour 
-æ, (0+), -æ is very useful, e.g. asymptotic expansions in 
cable problems. 

Bromwich followed up his first paper in 1916 by a number of 
others dealing with various problems in mathematical physics. 
These contributions to the subject appeared from time to time 
during the years 1919-30. 


* See McLachlan, Phil. Mag. 25, 261, 1938. 
t Titchmarsh, Introduction to Theory of Fourier Integrals (1937). 
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Heamside’s procedure and the Mellin inversion theorem. 


Those who have used Heaviside’s operational procedure are aware 
of the possibility of pitfalls (Bromwich called them slippery places). 
Books have been written for technicians by technicians in which this 
procedure is crystallised, so to speak, into rules, and in a number of 
practical applications they are easy to apply. But the blind usage 
of a technique (which is never justified by formal processes and 
logical argument in the said books) invariably stimulates a psycho- 
logical reaction, and the technician often finds himself in a dilemma. 
A number of engineers have informed me that one can never be 
certain of the Heaviside method, but it seems to be satisfactory 
provided the result can be checked in some other way! There 
are, however, problems whose solution by expansion in series and 
“ algebrisation ” in the Heaviside manner are unsatisfactory, the 
form of solution throwing little or no light upon the physics of the 
situation.* 

Attempts have been made hitherto to justify the Heaviside pro- 
cedure formally, without recourse to complex integration, and to 
show that if p=d/dt, then p”, where v is real, can be interpreted in 
terms of the gamma function. But such theses do not carry con- 
viction. 

By using the Mellin inversion theorem, as in the example given 
above, the operational solution of a problem is found as a function 
of p, i.e. d( p), the former being a mere number whose real part > 0. 
There is no question of p being identified with d/dt. Having got 
the operational solution, it is interpreted in terms of t by aid of 
integral (11). It happens, however, that the same operational solu- 
tion is obtained by writing p for d/dt, and the same function of t 
appears when the operational solution is algebrised in the Heaviside 
manner. Interpretation by the contour integral frequently gives 
the result in the form of a tabulated function, whereas algebrising 
may yield an unrecognisable infinite series. 

We may say that Heaviside’s operational procedure “ happened ”, 
as do many great inventions and events in life. They are none the 
worse for this. If experimental results could always be deduced 
a priori, research would be unnecessary, and this would deprive the 
scientist not only of his joy in life, but of his livelihood! The position 
respecting operational procedure may be summarised by saying that 
“ something ‘ happened ’ to Heaviside which made Bromwich and 
Wagner ‘ function ’,” to the great advantage of mathematicians and 
technicians. With the march of time, it is necessary to replace old 
methods by new in every department of scientific activity, wherever 
this is feasible. Moreover, the time is now ripe to remember what 


* See McLachlan, Math. Gaz. loc. cit. 


t As an example of this we can take the case of 1/(1 +ap2). The contour 
integral method yields 1— et/% (1 — erf V't/a?), but the series for this would not be 
easily recognised. This occurs in a problem on elasto-viscosity of materials like 
flour-dough or gelatine solution. 
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the illustrious Heaviside did, but to forget how he did it. So far 
as our present knowledge goes, the best procedure seems to be the 
application of the Mellin inversion theorem, as in the example given 
above. * When this is used, “ p ” is merely an inversion or trans- 
form parameter. Accordingly the nomenclature “ operator ” and 
“ operational ” is inaccurate. But in this matter we have to con- 
tend with usage extending over half a century, and an attempt to 
rechristen using, say, ‘ Parametric °” or “ Transform ”’ Calculus, 
would doubtless be greeted with ironical cheers. In any case a 
change of nomenclature would be unwise before being agreed upon 
at an International Conference. At some future date, when the 
Mellin theorem holds the operational field, the “ young idea ” will 
wonder how the word “ operational ? came into being to describe 
something entirely different. Some of the present generation are at 
a loss to explain why Thomson and Tait (T and T’) perpetrated the 
ill-founded title “‘ Spherical harmonics ”’. Musicians and acousticians 
should be sole concessionaires of the word harmonics ! 

Finally, we shall compare operational procedure with a time- 
honoured engineering process. Suppose we have a large ingot of 
steel and want to make an elaborate framework therefrom. By 
using hammers, chisels, hacksaws, drills and other engineering 
appliances, the ingot can be shaped ultimately into conformity with 
our design. This may be regarded as the sledge-hammer or battle- 
axe method. Ifthe ingot is melted and run into a sand mould (made 
by aid of a wooden pattern of the framework), the required result is 
obtained very quickly and much more simply than by the sledge- 
hammer way. We may conceive the moulding procedure to be akin 
to the operational method in mathematics, the sledge-hammer 
method comparing with analysis prior to the operational era. 

Transformation is made from the solid in a certain form to the 
liquid state, which after being given a predetermined shape is 
allowed to regain the solid state. The dependent variable y in a 
differential equation is expressed in terms of a parameter p, the 
resulting expression is interpreted, and the dependent variable y is 
obtained explicitly in terms of the independent variable t. 

N. W. M'L. 

1189. To offer the reader that sentence (from an Act of Parliament)... 
would be as little reasonable as the revenge offered to a village schoolmaster 
who, having complained that the whole of his little treatise on the Differential 
Calculus was printed bodily in one of the earlier editions of the Encyclopaedia 
Britannica, was told that he was welcome, in his turn, to incorporate the 
Encyclopaedia Britannica in the next edition of his little treatise.—J. Hill 
Burton, The Book Hunter, p. 144. [Per Mr. J. B. Bretherton.] 

1190. Between these two parallel lines a gulf still opens. It is far narrower 
than it was ; but still as profound and still unbridged.—Julian Huxley, Essays 
in Popular Science (Pelican Books), p. 138. [Per Dr. C. Fox.] 

* The most powerful method for solving technical problems of various kinds is 


that of integral equations. The Mellin theorem provides a ready method of solution 
for a certain class of problem. 


